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ABSTRACE 


Effects of a turbulent gas flow on the Stabile, gene 
glow discharge are investigated. While it is known that 
turbulence affects the stability of glow discharges, the 
mechanism is not clear. Primarily, the sropien lies in the 
very long characteristic times of turbulence phenomena when 
compared to the glow discharge instability times. 

A mathematical model is developed which is solves numer- 
ically for the ambipolar diffusion and an overall cheng, 
equation in an unsteady, cylindrical coordinate systen. 
Strong perturbations of the electric field are introduced 
which disturb an otherwise stable configuration and the 
effects of turbulence on the time evolution of the pereueeae 
tion are observed. 

It is shown that the modification of the ambipolar 
diffusion coefficient and the thermal conductivity is a 
reasonably sufficient model to introduce the turbulence 
effects. It is found that the charged particle density is 
effectively the sole source of heating in the gaseous environ- 
ment of the discharge. It is then shown that turbulence 
acts to suppress the temperature instability introduced by 


the discharge streamer. 
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Poe rvcne our TON 


BencegedsssucmM eas Nitrogen Of argon, Or a mixture of 
Peiseomemech asewadlt, 1S sWbjected stouvan electrical potential, 
a sequence of phenomena take place as a result of a continu- 
ously increasing applied voltage and/or current. Initially, 
for a small voltage, a very weak current on the order of 0.1 
picoamperes develops. This current will increase slightly 
with increasing voltage until all the naturally occurring 
charge carriers are being utilized. A current saturation 
condition now prevails until the applied voltage reaches 
about 6 kV in air at which time a further voltage rise 
rapidly increases the current flow to about one picoampere. 
This is referred to as the Townsend Discharge Region and it 
is due to new charge carriers being created by collisions 


between the highly energetic electrons (fast moving) and the 


neutral (slow moving) atoms. This process can have several 
effects on the individual gas molecules. Essentially, it 
can excite them, dissociate them, and/or ionize them. This 


ionization effect is higher for lower gas pressures since 
the farther the electron travels between collisions the 
higher will be its kinetic energy which is imparted by the 
acceleration produced by the applied electric field. When 
the electron's energy is greater than that required for 


10onization of the neutral gas particles the collisions will 


el 


produce additional electrons which in turn gain energy, 
produce collisions, and therefore, increase the total number 
of electrons present. This process is called an electron 
avalanemie. 

From this point, slight increases in voltage produce 
Current increases of several orders of magnitude at which 
time the current becomes self sustaining. Reducing the 
external circuit resistance now increases the current flow 
and a steady glow is produced in the gas. This glow- 
discharge is the familiar phenomenon which occurs in neon 
lights and has other practical applications suche ase 
pumping scheme for gas lasers. 

Finally, when the current is increased still further, a 
highly ionized channel is formed which collapses the glow. 
If there is sufficient current carrying capacity au eea— 
external circuit a steady current is established which 
represents a plasma state called an arc. 

A full understanding of the arcing mechanism is diffi- 
cult due to the extremely short formation times involved. 
However, Since delaying the point of transition to arcing is 
a highly practical problem, a method for modeling 14) 2 
needed. A blend of theory and experiment is usually re- 
quired. Such a model could then be used to Study seme 
effects of the arc on and by Jes envVinenmenc. 

It has been observed that this glow-to-arce transition 


can be delayed by utilizing gas dynamic stabilization 


(ky 


techniques in flowing gas systems [Ref. 1]. The purpose of 
this work is to establish the extent of fluid dynamic 
stabilization that results from the introduction of varying 
levels of fluid dynamic turbulence into the flow. It will 

be shown that this method has beneficial effects on the 
Stability of the glow discharge process and on the temperature 
distribution in the gas produced by the presence of discharge 
streamers. For further discussion and background information 
merating £O 10nization and the breakdown processes refer to 
any standard reference in the field or to Wallace, R.J. 

PReft. 2]. 

An introductory skill level in atomic physics plus 
fundamental proficiencies in mechanics, electrical princi- 
ples and numerical analysis should provide the reader with 
adequate preparation for understanding this material. 

Ti weoingelpbes discussed here have many practical engi- 
neering applications. For maximum power output of an elec- 
trically pumped gas laser, such as a carbon dioxide laser, 
the glow discharge must be held essentially at the limit of 
arcing. Since turbulent flows tend to delay this glow-to-arc 
transition, the laser can be allowed to operate at either 
the highest possible power level or with a greater margin of 
safety at lower power levels. Another application where a 
stabilized high frequency electrical spark discharge could 
be used as a means for injecting heat into a propulsive duct 


(like a pulsejet) has been suggested by C.E. Tharratt [Ref. 3]. 


i 


This type of a device is theoretically capable of unlimited 
thrust. It would also have a low specific fuel Consumpbeien 
and could operate outside the audible range. The possibile, 
of nearly silent operation would be highly advantageous as 

an aircraft propulsion device for heme cure- 

Two major cases will be discussed. The first case 
concerns the behavior of a highly transient discharge in 4 
non-ionized gaseous environment such as the lightning dis- 
charge. This first case concerns a pre-glow discharge type 
phenomenon since the streamer propagates in a non-ionized 
medium. However, it serves to establish the effects of 
turbulence in a simpler situation in order to form a starting 
place for this research. The second case will look at dis- 
charge streamers ina partially ionized gas which might be 
used for the pumping mechanism of a gas laser. The method 
of approach to these problems will be developed in the next 
chapter where a mathematical model will be presented. The 
governing partial differential equations for the model pre- 
sented will be solved by numerical integration techniques 
utilizing a computer. A working Fortran program 1S presemead 
which can be used to investigate a wide variety of the aspects 
of this problem. Significant results are then presented in 
graphical format to demonstrate the stabilizing features of 


Eheerurbiwilent yellows, 
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II. ANALYSIS 


eee bE CONTINUITY EQUATION 

Highly ionized channels which develop at the final 
stages of a glow discharge just prior to arc breakdown are 
called streamers. The transition to arcing is then often 
due to streamer-breakdown. It should be noted here that the 
Ilyper OL Streamers to be analyzed are the non-self-sustaining 
discharge streamers that precede the arc, rather than the 
electron avalanche (onset streamers) type that occur earlier 
in the glow discharge phenomena. 

Streamers are seen to issue from either electrode and 
appear always to precede the breakdown. These streamers 
propagate virtually instantaneously for our purposes as an 
Mbeiatast LONiZation wave. They consist of a partially 
ionized plasma column for which the conductivity is dependent 
directly on the charged particle density within the streamer 
[Ref. 4]. It will be shown that the collisions made by the 
Sm@argeda particles within these streamers are the main soieee 
of heating of the surrounding gas. The value of this 
charged particle density may be determined by the applica- 


tion of the following continuity equation to the breakdown 


process. 
D 
A Zz 
Bee eitgy) — 0 se - ajn* - an” pea 


eS 


Equation (2.1) gives the time rate of change of Ehe "enaraea 
particle density, n, as a function Of radial post le. ean 

and time, t, for the decaying streamer. It is written in 
cylindrical coordinates due to the natural form of the 
streamer. A schematic drawing of the streamer geometry is 
depicted in Figure 2.1. The right hand sidowy aee 
Equation (2.1) consists of three classes of terms. The first 
term represents the effect of diffusion due to the niga 
density gradient within the plasma channel. De 1s the ambi- 
polar diffusion coefficient. This ambipolar diffusion resem 
in the electrons and ions diffusing out of the original 
channel at essentially the same rate. Briefly, the lighter 
electrons having higher thermal velocities leave the plasma 
ahead of the slower, heavier ions. This produces a local 
positive space charge which retards the loss of electrons and 
accelerates the loss of. the ions. The net effect is that the 
ambipolar diffusion coefficient is approximately twice the 
diffusion coefficient of the ions alone [Ref. 4]. 

Next consider terms 3 and 4 of the RHS of Equation (273) 
which represent two and three-body recombination respectively. 
Recombination is a de-ionization process which is then 
responsible here for the reduction in conductivity of the 
partially ionized gas. These processes occur when an electron 
and an ion collide at relatively low velocity and recombine 
to form a neutral atom. In order to conserve linear momen- 


tum a third body must be present. In two-body recombination 
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Figure 2.1. Streamer Geometry 


1k 


this third body is a photon and, therefore, this proceso 
sometimes called radiative recombination. For niu eos eee 
recombination the third body is usually an ion or neutral 
atom [Ref. 4]. Due to the large complexity of electron meee 
mechanisms accurate values for the recombination coefficients 
are difficult to obtain. The values used in thusmeme see 

will be discussed in the next subsection. 

The second term on the RHS is the convection term. This 
term is neglected since the characteristic time of convection 
(length/flow velocity) is much less than the characteris 
recombination times associated with terms 3 and 4 [Ref. 2]. 

While the convection term in Equation (2.1) can be 
neglected, the diffusion term cannot since the diffusion 
coefficient D, can increase by about three orders of magni- 
tude between laminar and turbulent flow [Ref. 2]. This will 
cause the characteristic diffusion time to be on the same 
order of magnitude as the recombination time and some rela- 
tionship will exist between diffusion and recombination which 
will affect the charged particle density distribution yea 
in space and time. 

Once Equation (2.1) is solved for n(r,t) the resulting 
distribution profile can be integrated to determine a 
conductance of the decaying streamer as a function of the 
streamer radius and time. This can be accomplished under 
varying turbulence level conditions by changing the erred 


tive value of the @QiffEusion Coct frelon 
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Equation (2.1) describes the lightning-discharge type 
phenomenon in a non-ionized medium. A more interesting 
situation is the case of the streamer as it appears in an 


1onized-sustainer discharge for the pumping of a gas laser 


Meet. 5]. For this case, Equation (2.1) may be written as, 
on ea 3 on 2 3 
a SS ( ray = vin - ajn” - aan (22.2) 


where now a new first term is added to the RHS which repre- 
sents a production of electrons and ions. The ionization 
coefficient, Ves is a function of the sustainer E/N and can 
be determined from Figure 2.2. The electric field referred 
to here is an existing axial electric field, that is, the 
Boemping field” that maintains the glow. Equations (2.1) 
pum 2.2) are non-linear partial differential equations 
which will be solved numerically by the use of a finite 


difference approximation scheme. 


Eve LHE ENERGY EQUATION 

In addition Ecco hiw ity equation, a Suitable 
energy equation must be considered in order to describe the 
Bemperature distribution introduced into the gas by the glow 
mesecharge streamer. For thermal stability to exist, a 
balance must be created such that the time rate of change of 
the internal energy is zero because the loss of heat through 


thermal conductivity and the addition of energy from Ohmic 


159 
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heating must balance each other. The following equation 
represents this relationship. 
De = xa = 
ai ° = . De 
P De E 5 V gq (225) 
The first term on the RHS of Equation (2.3) represents 
Joule or Ohmic heating where, 2) is the electric field 
strength and, 5, is the current density. The second term 


is the heat flux loss due to thermal conductivity and 


radiation. The total derivative, 
De Je de de de 
a 7D | cna ip eee Je ore ; 
Die 0 a IS yaoi : 6 08 7 Zand Z (2.4) 
reduces to, =, since the model under consideration here 


feeovs £Or Conduction only in the "z" direction. It is 
further assumed that the change of the internal energy along 
the length of the streamer is small, the radiation emitted 
primarily perpendicular to the streamer centerline has a 
negligible effect and, due to the short characteristic times 
involved, the heating occurs at an essentially constant 


volume. With these assumptions in mind, Equation (2.3) 


becomes, 
oT 7 2 
OC ae KV T = one (2S) 


miaere, Pp, is the gas density, o, is the electrical conduc- 


tivity of the plasma (recalling that j = cE.) and K is the 


Za 


thermal conductivity. Equation (2.5) gives the temperataeer 
T, as a function of the radius andecmee 


The conductivity, o, 2S gGilvoomers 


o = eniy (2 Jon 


where, e, is the elementary electronic charge in coulombs, 
nes 1s the number of electrons and, Wes is the electron 
mobility. The mobility is also a function of E/N and is 


defined as the ratio of the electron drift velocity, Le 


Vat 
the applied external electric field. The value of V5 is 
given in Figure 2.2%as a function oOfeE N- 

It will be shown later that for the OES term from 
EQuUatrOn | 2.5), ee is essentially constant and since o is 
directly proportional to n the solution of “Equation (2s 
represents the single source of heating that drives the 
temperature distribution opened from Equation (2.5) 2 eae 
consequence of this is that a mathematical model describamg 
the temporal and spatial characteristics of the streamer can 
now be developed. Furthermore, the effects of turbulence as 
a stabilizing medium can then be studied since the thermal 
conductivity, K, and the ambipolar diffusion coerriciemee 


D are the only parameters affected by turbulence in Equations 


a 7 


(251) (252) eandet2.5)— 


C. SOLVING THE VEOUATIONsS 
An effective means of generalizing a problem such as 


this is to non-dimensionalize the equations. To begin then, 


22 


after some algebraic manipulation, Equation (2.1) can be 


merpeten in the form, 





& D Cia) 
S A2 3 a3 
—— 22 - (—45)4tr B= wn*® - = 2 (2.7) 
ZO he aes or ox Z 
© 2250 
where n= Ayn. -and r= ie Vem ct. are the initial 
Oo O Oo O 


values for the charged particle density and the streamer 
radius. The following additional definitions will prove 


useful in this model, 





pone 
SS 5 (Zero) 
a 
Clee hl 
Qa = 3.9 (29)) 
nD 
t = = akon 
R 


The parameter, 8, 1S an inverse measure of the turbu- 
lence level since the ambipolar diffusion coefficient is 
smaller in laminar flow (large 8) and larger in a turbulent 
flow (small 8). If one neglects the diffusion term in 
Bepeeron (2.1), then for two-body recombination, direct 
integration yields, 


re 


a O 
a ort (2.11) 
2 0 


t~eselon loss tormula [Ret. 2] then shows that the 


Characteristic recombination time, Tre is simply l/fa5n. - 


a3 


Incorporating the above definitions allows Equation (277) 


to be written in the non-dimensionalized form, 


on = i care oi, = ase = ane (2.14 


n~ 


st PP ar or 


which lends itself to a numerical eo wRei1en by substituting 
the appropriate finite difference approximation ferme 
(Ref. 6] for the differentials. After some manipulation, 


Equation (2.12) can be approximatedge. 


“a _ “~ I _ w~ 
Nee oaal = LAtIZ(A +B) Cj} + ee (2. ie 
where, 
A = nae Mitl,3 ie = eee (2, as 
oe ee . ae * ( 2:aes 
ne aera i-l,j 
ee a3 
Ca Ons + ee (2, lay 


In this scheme, ‘'i' represents position perpendicular to 


'j' represents relative position 


the streamer centerline and 
in time after some initial condition. Equations Za 
through (2.16) form the basic algorithm used by the Fortran 


program discussed in Appendix A. 
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The accuracy of the charged particle density obtained by 
this type of differencing scheme depends on the step sizes 
At and Ar. The steps must be smaller than the smallest 
characteristic time of the phenomena described. Therefore, 
At is taken as 0.1 or in real time 1,/10 [Ref. 6]. Primarily 
for programming convenience, Ar is also taken as 0.1. Addi- 
tionally, the value of g is limited by the following stability 


Semcerion [Ref. 2]. 





gp > so Gay 
pee 
Returning to Equation (2.2), it is seen that the addition 
of the ionization term to the RHS requires a modification 


to Equation (2.16) which now becomes, 


Cc = c + aia -on. .] (er) 
rd Lo a 


Equation (2.5) can be non-dimensionalized in a similar 


manner by writing it in the forn, 








A 2 
L ro ee Ob 
a 
(Sehr SH - (2.19) 
A Jt, pC ry ioe FC ox V 


where TA is defined as, 


Me (2.0) 
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A check of Equation (2.20) reveals that Ta has the 
units of seconds and, therefore, is a characteristic time 
for the energy equation. If one now defines the following 


additional relations, 


ct > 
(Wt 


(2.25) 
A Ty : 
a e lt 
| = rT (2,225 
O 
OE r¢ 
O 
B =. Jace (2.239 
A KT ' 


then, it can be shown that, Equation (2.19) becomes the 


following non-dimensional relation, 








Tr iL a 
ee Sr SE) SF) (2.24) 
Jt, Bar ox ox Oo Prag 


where, ay is the initial temperature, which for this work, 
is taken as 300° Kelvin. It should be clear here that, Bae 
1S now a non-dimensional parameter whose value depends only 
on the thermal conductivity, K. This parameter then serves 
the same function for the energy equation as does 8 for the 
continuity equation, that is, it allows the influence of 
turbulence to be introduced into the model. The value of 
the electric field strength, Bo from the RHS of Equaeirer 


(2.24) can be estimated by assuming that the initial value 


of the externally applied field, Bo! is given by [Ref. 8], 


20 


E = E+ — (32 225) 


A reasonable value for 2/2 might be on the order of 


meooe Q/m, for which Equatiom (2.25) becomes, 


= I 
ES =e (lee 7 (10,000) } (2ae7 6) 


where I/E for the streamer can be shown to be, 


e(u_N) x 
CT _ e 
aaa (2a) ) 
0 
enc! , 
Ee 
N = oF 22 e) 


"N" is the total number of particles present in the gas at a 
pressure, P, and a temperature, T. The Boltzman constant is 
mae cr= 25 7I/ K. The integral in Equation (2.2/7) is a 

measure of the conductance of the streamer and is evaluated 


in Chapter III by integrating the charged particle density 


PieeriouLiOn Obtained from Equation (2.1). The initial 
value of this integral is simply more. Using estimates 
3 


of n_ and ry to be 3.25E18 m~ and 4.0E-4 m respectively 


O 


[Ref. 7], Equation (2.26) can be solved for E. becoming, 


2 


— =. ee at 
” T+ 0.0000 mae (2.29) 


where the value of HN in Equation (2.27) can be eSeaumec 
From reference tables to be 5.257E24 (v-m-s)~+ for nitrogen 
Results to be discussed in Chapter III show that the value 
of the Taeeceml drops rapidly to a small and nearly steady 
value (on the order of 0.01) such that the denominator of 
Equation (2.29) approaches unity in very sheresemecr 

The results of the previous paragraph then show that the 
RHS of Equation (2224) is \sampoky n Since Ox is (en u/enou,) 
and Bie 1S approximately one. So now finally, Equation 


(2.24) can be written as, 








—] = on (2.308 


Equation (2.30) then shows that the temperature distri- 
bution due to the presence of the streamer is affected only 
by the charged particle density. 

As for the continuity equations, Equation (2-30) "eam 
solved by introducing the appropriate finite difference 
formula as follows (after some rearranging of terms), 


n~ 


See 
ee Oe ede + eee + I ae (2.308) 


J > 


where, 
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_ il 
and, 
Se ite (ees) 
eS a+1,) i-1,3 : 


The Fortran program discussed in Appendix A solves 
meameeton (2.31) simultaneously with Equation (2.1) or (2.2) 
depending on the case being considered by using the new 
value of n calculated from each step as an input to the 
calculation of the new temperature value. The time step for 
the temperature calculation is adjusted such that the real 


time step for the two equations is equal. 


D. MODELING THE STREAMER 

Mae results of Chapter Ili show that appropriate selec- 
tion of the controlling parameters for Equations (2.2) and 
(2.30) will produce a stable system in space and time. Then 
the question is: what effect does the momentary appearance 
Sia Streamer have on this stable situation? If a local 
Perturbation of the electric field occurs, then the ioniza- 
tion coefficient will momentarily be raised which will in 
turn raise the charged particle density in the region of the 
perturbation. Figure 2.2 provides the data for estimating 
memmagnitude of this perturbation by noting the relative 


increase in E/N for a corresponding increase in v,/N. The 


Oe) 


ratio of these values can then be used to determine a 


Simple multiplier (labelled Z and zy) for the n terms 


appearing on the RHS of Equations (2.52) and (2730) com oa 








late the effect of the streamer. Therefore, one can write, 
on - ee ae es = i = a - oe (2.349 
at B,ror ei 
and 
oe _ ae ee ae a Zn (2,30 
ot, Bar au 
where Z and aN are the multipliers discussed. For wen 


finite difference approximations, these two additions are 
easily handled for Equations (2.2) and (2.30) respective 


as follows. 


C. = ine rar ee (2.363 
ma = m™N 7. m™ mas 
eet = 7c ee eae i ee + Tas (2.398 


The Fortran program incorporates this concept by utilizing 
a separate subroutine called "STREMR" which is identical to 
the subroutine "DENSTY" except for the addition of Equations 
(2.36) and (2.37). In this scheme then, the streamer Gs 
modeled by calling up this subroutine when desired and then 


observing the reaction of the system over time. 
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A. NON-IONIZED ENVIRONMENT CASE 

In this section the effects of turbulence on the 
streamer in a non-ionized environment are considered. 
meamaeton (2.12) applies to this situation for which Figure 
3.1 shows a typical plot of the charged particle density 
versus radial position perpendicular to the centerline ina 
fully laminar flow. Note that the abscissa and the ordinate 
are expressed in the non-dimensional coordinates of r = r/X 
and 7 = n/n. respectively. There are five curves plotted 
which show the charged particle density distribution profile 
from the initial conditions at time zero in steps as specified 
in the legend to the profile at 3000 recombination times (TR). 
"TR" is the computer's version of "Tp and will be used here 
in the text. Now, for this plot, and except for Figures 3.11 
mmc. t2, £Lor all the figures in this section, no = 3 2 ees 
ine and a5 = Zao B= 12 W/E [Ref. 7]. These values then 
equate to a recombination time of 0.16 microsecond, and so, 
1OOOTR is equivalent to 0.16 millisecond in real time. 

In the laminar case depicted in Figure 3.1 it is clearly 
Peon that the centerline charged particle density drops off 


Mieelme In accordance with Equation (2.11) which in non- 


dimensional form would be, 
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For example, at l1000TR, Ne = 1.0E-3. This decay is due 
Permarily Eoucecombination only since, in this case, 

beta = 50,000, and therefore, the diffusion term becomes very 
small. For longer times, diffusion comes into play and some 
spreading of the charge profile is seen. 

Figure 3.2 shows the effect of introducing increasing 
amounts of turbulence into the flow. This plot compares the 
different charged particle profiles as they would appear at 
Mes2 msec £Or beta values of 50 (fully turbulent), 500, 5000 
and 50,000 (fully laminar). Note that the 2000TR curve from 
Figure 3.1 is identical to the beta = 50,000 curve in Figure 
Bere. The effect of the turbulence is clearly seen to spread 
out the profile and to lower : on the centerline as well 
[Ref. 2]. 

As was mentioned previously, the conductance of the 
decaying streamer can be determined by integrating the 
Sem@aeged particle density over the radius. This can be 
derived by "inverting" Ohm's law and introducing the conduc- 
tivity expression for a partially ionized plasma for which 
the current can now be expressed as [Ref. 4], 


Ben A 
2 


IT = GV = (—)V = ( )V (3a) 


mov. 
ee 


where G is the conductance, or the reciprocal of the resis- 
mance, in the familiar relation, V = IR. 'o' is the 


conductivity, 'A' the cross-sectional area and '2' is the 
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length of the streamer (or the distance between the elec- 
meeades im Figure 2.1 for this case). Additionally, m and 
Vv, are the electron mass and the electron collision frequency, 
respectively. 

For the decaying streamer, the radius, and therefore the 
area, is changing but the length is constant. Since all the 


other factors within the brackets of Equation (3.2) are also 


constant, 'G' can then be expressed as, 
ee Senn eR 
G = tm vy i eee Tee Guin G3) 
ee 0 


where gamma is now defined as the value of the integral in 
the above equation. Gamma has the units of reciprocal 
meters, and is determined in the program, "CHARGE" (Appendix 
Peeeoy the subroutine, "CNDUCT." Gamma is then directly 
proportional to the conductance. 

Figure 3.3 shows gamma plotted over time out to LOOTR. 
Note that gamma drops in value by 90% in only about 10TR 
(1.6 microsec). After this time the value levels off to 
pout 0.04 at LOOTR and then beyond the range of Figure 3.3 
meer value drops very slowly to about 0.001 at 10,000TR (1.6 
msec). 

pegume 3.3 1S for fully laminar flow. Introducing 
turbulence into the flow causes a very slight increase in 


gamma over all time for the complete streamer. Figure 3.4 
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shows this effect at lOOOTR. Figure 3.5 shows the value of 
gamma within the original streamer tube channel only. Here 
it is seen that for small beta (high turbulence) the conduc- 
tance is lower than in the laminar flow condition. 

Referring back to Figure 3.2, and with these conductance 
characteristics in mind, it seems rather apparent that one 
"effect of the turbulence is to delay the possibilicy sei 
transition to arcing since the conductance is lower along 
the original channel. Another consideration is what happens 
to a follow-on streamer appearing after the charged particle 
density profile of the first streamer has died away for say, 
2Z000TR as in Figure 3.2. Ina laminar flow there would be 
a distinct path of least resistance along the original 
Channel. However, in the turbulent flow, the profile is 
Spread out over a much larger radius and then, this "“follow- 
on streamer" essentially has no favored "route" for arcing. 

Equation (2.30) showed that the charged particle density 
was the sole source of heating in this model. Figure 3.6 
shows the result of solving Equation (2.30) ina fully 
laminar flow for an E/N Of 22 63-ba V-cm", at one atmosphere, 
and 300°K for an electric field strength of 4.0E4 V/m. Note 
that the temperature rise amounts to less than about 0.05% 
in 1000TR and then rises very slightly over time eventually 
peaking out at about 4000TR and then falling back slowly. 
This would seem to be a logical behavior for this case since 


the charged particle density that is heating the gas is also 
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dying away. The question is then; what happens if one 
increases the electric field? 

Paice se ey weineutiecew 0 eslliowetme effects of increasing 
Ey from 1.0E6 V/m (a typical value of interest for lasers) 
fono.0E6 V/m. Note that Pn changes for each of these cases. 
mimes aS due to the fact that the electron drift velocity is 
Pmetetion Of E/N, which is plotted in Figure 2.2. Since the 
drift velocity changes, the electron mobility changes which, 
in turn, produces a new conductivity value for each different 
electric field strength case. Recalling Equations (2.20) 
and (2.23), one can see that each case has a different 
characteristic time and a different Bae For increasing Boe 


tT. gets smaller while Bn gets larger. 


A 
Riemmali Gdietficulty with this 1S that the time step for 
the energy equation solution must be adjusted to ensure that 
miewconcinulty equation is “feeding” the charged particle 
density values into the energy equation at the correct rate. 
For the continuity equation, the non-dimensionalized time 
Seep is 0.1 (T/10 in real time). Therefore, the true time 
for one time step (one "DT" in the program) is simply 
Piet. Similarly, the true time for the energy equation is 


R 
POTA Ne It follows then that the value needed for DTA is 
simply (DT XTp)/T,- 

PaeUibadgmre Pigures 3.7 through 3.10, it can be clearly 
seen that increasing ae increases the temperature rise, as 


would be expected. The effects of longer times is the same 


femtme Case for Pigure 3.6 but requires much longer computing 
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times to reach the peak and to begin tfalVing bach 
effect of turbulence on the Cases Of interest, say 40084 

V/m to 1.0E6 V/m, is to essentially eliminate the temperature 
mise. 

One remaining question for this non-ionized case is what 
effect would increasing ny have? Figure 3.11 shows that at 
ES = 1.0E6 V/m and Ls increased by a factor of 10 to 3.125n1 97m 
produces a more squared off profile than that depicted in 
Figure 3.7 suggesting that the temperature dies rapidly 
outside of the immediate influence of the streamer. Figure 
3.12 shows the effect of another 10-fold increase in N. 

Here the resulting temperature rise drops off at an even 
steeper rate at 2000TR than for Figures 3.7 and 3.1ll. 
Recalling that TL = T/asno. it is clear that increasing ne 


R 


decreases the characteristic recombination time significantly 
which, in turn, lowers n such that there is a noticeable 
reduction in the temperature profile outside of this immedi- 


ate streamer influence. 


isa IONIZED ENVIRONMENT CASE 

In this section the effects of a glow discharge in an 
environment where a steady glow exists is considered. For 
this situation, Equation (3.4) is the applicable non-dimeneaenes 


form of the governing continudey equacien: 
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The temperature effects are then determined by Equation 

(2.30) where n tomOocained by Solvang Equation (3.4). The 
Earst question in this case might be; what effect does 
turbulence have on the temperature profile due to this 

steady glow discharge? At sub-atmospheric pressure, it is 
known that the charged particle density profile takes ona 
Bessel function distribution [Ref. 9]. Therefore, to inves- 
tigate this situation, Equation (2.30) can be solved inde- 
pendently from Equation (3.4) by calling up the subroutine 
BeENSTY” ain the program and sending it a constant, for all 
time, value for — The Bessel function serves as the initial 
condition profile for the value of 7 ih Seis teen them, 
labeled "Nl" in the program, allows : to go to zero at r = 4.8. 
Figure 3.13 shows the results of such a case for strictly 
memenar Flow. The values of the controlling parameters for 
megures 3.13 through 3.16 are: n= Se 258 18 ae 
a. = eon 12 m>/sec and EG = 1.0E6 V/m. It is clearly evident 
that the temperature for such a laminar case would rise to 
approximately 5400°K (recall that T/T = lis eguivalent to 
300°K) in only about 1000TR. Since molecular nitrogen (N.) 
begins to dissociate at about 4000°K, the situation depicted 
in Figure 3.13 is undesirable. However, introducing turbu- 
lence into the flow tends to significantly reduce this 
temperature distribution. Figure 3.14 shows that a moderate 


degree of turbulence can be selected which, essentially, 


holds the temperature profile near the laminar profile without 
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going unstable as in Pigure 3.13. Pigqures 3.15 a0 ee 
show the effects of increasing the turbulence level still 
Further. Note that in Figure 3.15, the temperature promi 
is nearly "frozen out" at about 2400°K, while Figure 3.16 
is holding steady ae around 220071 

This asymptotic behavior of the temperature rise in 
reaching this stable profile appears to be related to the 
Characteristic time associated with Equation (2.30). Observe 
that for all three figures (3.14, 3.15 and 3.16) the tempera= 
ture rise stabilizes between about 2000TR and 3000TR. This 
effect will continue to be evident in the remaining cases to 
be discussed. The parameters for the case presently being 
considered are such that the true time of 2500TR is 0.4 msec. 
This is approximately 4.5 times the characteristic time step 
in seconds for the energy equation, which is consistent with 
the concept of a rise time in an asymptotic system. 

The final case to be discussed here is perhaps the most 
interesting. This will be the case where there is a steady 
level of ionization in the gas and a perturbation of the 


electric field is introduced to model the eGrieccts oreaan— 


streamer. For the steady state situation at the centerline, 
3 = O, sO one can write Equation...) 
D 
aoa on _ = _ 2 3 
ze ap \F ar ¢ = 0 Vas Soe: a3n¢ (3.359) 


Neglecting the three-body recombination term (a. is on the 


Order tor 2). 0n—34 m°/sec), and Oe iNg@imm=ie a on the centerline 
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momalseeZero, tlie 1OniZation term is then equal to ee 
Since Figure 2.2 gives V./N, the charged particle density 


for the steady state condition can be obtained by, 


Ve vy 


ee 2 1) (Nv 
et = a5 = ( Nn? ‘a5? (3 a6P 
where N is given by Equation (2.28). Figure 3.17 shows this 


Steady state charged particle distribution where, for this 
model, the glow discharge is contained ina region of a 
radius equal to 10 ro: At Lo! the boundary condition is set 
eweh that n is zero at the "wall." The left hand boundary 
condition assumes that the slope of the charged particle 
density is zero. See Appendix A for details concerning how 
the program meets these boundary conditions. Note that in 
Figure 3.17, the steady Bee eeeon ancien is reached in less 
than LOOTR, and is constant at the longer times. The first 
vertical line, at r/X = sie) ms the initial condition 
distribution which assumes = Oe eUlinat© r = 1.0 and zero 
beyond this point. 


Mtewsteady State solution for Equation (2.30) is given 


Mmmequation (3.7). 


“a 2 2 2 “A 


ot i nb 
T = (8.n) re te By aan + V5 (37) 
Geen boundary conditions are T = Bo ese e—o) “and ot =a) ac 


r= 0. Figure 3.18 shows this result for ES = 1.47E5 V/m. 


For this case By = 0.2219, and so at the centerline r = 0 
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and substituting in the appropriate values into Equation 
(3.7) One ebtarns T = 6.548 which is confirmed at r = 0a 
Figure 3.18. Note that this steady stave condi elena 
reached at approximately 2500TR which is about the same time 
period required for stabilization from the previous case 

for the Bessel function initial condition protile (Prague. 
3.14, 3.15 and 3.16) and for the case showing the temperature 
profile of a lone streamer in a non-ionized environment 
(Frqures 3.6 throughwea. 2) 

Figures 3.19 through 3.28 show the results of introdueaiag 
the electric field perturbation into the model. Each ploer 
shows the situation of the streamer appearing in a steady 
state turbulent environment and lasting for increasing 
periods of time beginning at 5 recombination times and being 
incremented by 5TR up to a final value of 50TR. In real time, 
this would be equivalent to 0.6 msec to 6.0 msec. For these 
plots; oe = 1.47E5 V/m which produces an E/N of 6 20ESay 
V-cm*. This case was chosen for example purposes since the 
steady state temperature profile is convenient at a centerline 
value of T7 ed = 6.5 and the rise due to the perturbation 
reaches a maximum of about 13.75 (4125°K) at 50TR. Each figure 
shows five curves, 3 are solid and 2 are dashed. The solid 
lines represent the initial temperature distribution (T = 1.0)9 
the steady state profile, and then, the profile of the streamer 


at the end of its stay time. The dashed lines are the decay— 


ing, or transient, profiles produced after the appearance ma. 
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the streamer. Note that the centerline temperature drops to 
about half of the rise value in only 20TR and has nearly 
reached steady state in another l100TR. 

Figures 3.29 and 3.30 show a case for a higher E/N value 
oe. OE-1/ one Here the steady state temperature profile 


has a centerline value of about 10.0 and then a streamer "on 
time" of only 8TR sends the profile above 4000°K. Cases for 
smaller E/N values produce lower steady state temperature 
profiles. The streamer effects are similar except that 

much longer "on times" are needed to produce much of a 
temperature rise. This appears consistent in the model 
Since at the lower temperature it would take a much larger 
perturbation to upset the steady state condition, while at 
the higher energy levels a comparatively smaller perturbation 


is all that is required. 
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This thesis presents a mathematical model for describing 
the effects of varying levels of turbulence on glevean 
charge phenomena. The development is concerned only with 
the aftereffects of the streamer, rather than its propag@as 
tion. The program that is presented in Appendix A is 
intended to allow the investigation of a large variety of 
cases of interest relating to the glow discharge. The cases 
discussed in the results chapter are those that appeared to 
develop a logical sequence of thought in attempting to 
answer the question of how turbulent gas flows can and do 
affect the stability of a glow discharge and the problem of 
Gelaying the glow-to-arc transition. 

The results clearly show that the modelled turbulent flow 
does indeed have strong beneficial effects on the stability 
of the system. Additionally, the mathematical model success- 
fully predicts within a reasonable level of uncertainty a 
region of stability that is consistent with the normal oper- 
ating ranges of, say, an electrically pumped gas laser 
[Ref. 10]. The model suggests that the region for which the 
system maintains thermal stability is critically dependent 
on E/N. The research showed that a steady state temperature 
increase from 1200°K to 3000°K results from an increase in 
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One important question still remains however. That is; 
is there a particular level of turbulence which would have 
the most beneficial effect on the streamer-initiated break- 
down problem? This work suggests that the best situation is 
Simply that the most desirable level is the highest level 
that can be maintained. However, the non-linear nature of 
the equations suggests that some optimal level may exist 
[Ref. 2]. Experimental work is probably the best means of 
answering this question. 

The results presented here by no means exhaust the aspects 
of this problem that could be investigated. The intent of 
the program is that it can be used as a tool to analyze any 
desired case. Some logical next steps might be to expand on 
the streamer modeling aspects to investigate what happens when 
Bmeeel iow-on streamer appears before the initial one has fully 
decayed back to the steady state condition. The "Z" and "2," 
arrays could be set up to have two or more streamers appear 
at any desired location. Also, the values of "Z" and "Zn" 
used -in this model are somewhat arbitrary and probably have 
to be specified for each individual E/N case according to 
experimental observation. More research into the actual 
magnitude of the electric field perturbation may allow the 
model to more accurately predict the areas of instability 
with streamer on times that are more in agreement with 
experiment [Ref. 5]. 

Mathematically, the analysis of this problem consists of 


solving the continuity equation for three cases. The first 


Te 


is for the RHS of Equation (3.4) equal to less thai =e 
This is the non-ionized environment case. The second is for 
the RHS equal to zero. This is the steady state case whies 
is what was applied to the ionized environment cases discussed. 
The third case, which was not applied in this thesis, Sipe 
the RHS of the continuity equation to be greater than zero. 
Investigation into this third case could possibly  elanime. 

some aspects of the problem. And finally, the energy equation 
may have eigenvalue solutions which might produce results of 


interest such as additional regions of stability. 
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APPENDIX A 


THE CHARGE PROGRAM 


The Fortran program "CHARGE" is intended to be a working 
tool for analyzing the glow discharge problem. It evolved 
from a previously written program which solved the non-linear 
continuity equation (2.1) [Ref. 2]. It is now arranged wane 
somewhat different form. The program consists of a main 
program section which calls up any needed combination of the 
seven subroutines which perform the desired tasks. As cur- 
rently written, it can compare four separate cases plus it 
will display the initial conditions of either the charged 
particle distribution or the temperature profile. Individual 
test cases are set up by commenting out (placing a "C" in 
column 1) those commands not desired. After the main program, 
there are two sections. The first provides commands to print 
out the data and the second provides commands to plot the 
results by utilizing the "DISSPLA" system developed by the 
Integrated Software Systems Corporation, San Diego, California. 

The program as it appears following this discussion is 
set up to produce the plot displayed in Figu@re 3.263 ea- 
main program begins by defining all arrays. FOr this ve7ane 
ple the arrays are set to 101 in the "i" direction which 
specifies the size of the array with respect to radial posi- 


tion. Since one DR, one radial step, is equal to ro/10, 
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pectaten UO es pecifies a distance of 10 E perpendicular 
from the centerline of the streamer. The "j" size is always 
kept at 2 due to the computing scheme used in the subroutines. 
After setting the array size, set the variable named 

"RADIUS" equal to the same value as the "i" coordinate above. 
In this case, RADIUS = 101. This determines how far in the 
radial direction that the differentiation is to procede. 

Next one must assign the desired values for BETA and BETAA 
depending on the case at hand. The BMIN, BMAX, BAMIN and 
BAMAX parameters apply when it is desired to transition the 
calculation from a laminar to a turbulent environment. This 
feature is designed to attempt to ensure that the calculation 
accounts for the fact that the streamer initially appears so 
quickly that the flow looks laminar to the streamer no matter 
what the level of turbulence is. Then, after a period of 
eeme, the charged particle distribution left behind by the 


streamer has existed for sufficient time to now "see" the 
turbulent effects. The time required for this transition is 
estimated by noting that the characteristic diffusion time 
for laminar flow is approximately 0.1 msec [Ref. 2]. Depend- 
ing on the case at hand, one DT is on the order of a few 
microseconds. Then, 100 DT (1OTR) is on the order of a few 
tenths of a millisecond. Therefore, for most cases, a 
transition of 100 DT is used. 


Returning to the example, one procedes by selecting the 


desired initial conditions. In this example, all of the 


ear 


"Nl" values are commented out since they apply to the Bessel 
function initial condition profile. Four series Ciysuprer 
tine calls are provided then to allow running Of v:eue 
separate cases. Each case assigns its respective results to 
the "VALX" (X = 0,1,2,3,4...) array £6, Sete rmigeee weer 

The initial conditions are assigned to VALO or VALTO# 
Then the first subroutine call will see the initial distribu- 
tion values of n while each successive call will see the 
results of the previous call. If it is desired to recover 
the initial conditions, call either “NZERO" or “TZERO @ac 
desired and the current values will then be reset to time 
zero. 

For this example now, "VAL1" is the result of two calls. 
The first is for "TRNSTY" to transition from laminavees 
turbulent conditions with time set to 101, then “DENSTY 
called and runs for 23,900 DT. The result is the steady 
state temperature distribution in turbulent flow at 2500TR. 
VAL2 is then the result of calling up the subroutine "STREMR" 
and running it for 500 DT with BETA and BETAA set to tie 
laminar values. VAL3 then results from calling up "TRNSTY" 
for 100 DT and "“DENSTY" for an additionailog@ eae 
turbulent values. Finally, VAL4 simply runs “DENSTY™ fea 
another 1000 DT in turbulent conditions to check the wdeaaw. 
ing temperature profile. Of course, any combination of these 
subroutine calls and time values can be selected to view the 


Situation as desired. 
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The results can then be printed out and/or plotted. To 
become familiar with the "DISSPLA" plot system, refer to 
DISSPLA Pocket Guide (current with version 9.0) by ISSCO 
eecphiics, 0505eserrento Valley Rd, San Diego, California. 

Some pointers of potential errors should be clarified. 
When selecting the desired heading lines, ensure that the 
last entry in the call to "HEADIN" correctly indicates the 
total number of heading lines. If no heading appears, this 
is probably the reason. Secondly, ensure that the third 
entry in the call to "CURVE" is equal to the "i" value on 
the arrays. This is the number of points (101 here) to be 
peoetca. And finally, the second entry in the call to 
"LEGEND" is the number of lines of legend needed, 5, in this 
example. 

wine Subroutines "DENSTY,"” “STREMR” and “TRNSTY” all 
contain essentially the same algorithm with modifications to 
perform their particular function. The basic idea here is 
ema the Calculation procedes from the first row of data 
and fills the second row with updated values. Then the first 
row of old data is replaced by the second row of new data by 
"SWITCH" and the loop increments one DT and procedes. With 
this scheme the program can be run as long as desired without 
using large values of memory storage. 

The boundary conditions for the problems are met by the 
routine. The left boundary condition assumes that the slope 


1s zero across the streamer centerline. So here, the new 
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an “* 


value at n(2,2) is assigned to n(1l,2) before the switch. 
The right boundary assumes that n is zero at the wall, that 
is station 101. This is met by limiting the maximum radius 
value to be one less than "RADIUS" (called "TEMRAD"). The 
routine procedes to "TEMRAD," then the last point seen by 
the calculation algorithm is zero. Similarly,,the tempera= 
ture of the right boundary is set to 1.0. 

To compute the conductance parameter, gamma, that would 
result from any particular a distribution, call "Gib 
after assigning the desired : distribution to some "VALX." 
The single value "NN" will be returned, which is gamma. 

It is hoped that the program is sufficiently self- 


explanatory. All variable, constant and subroutine names 


are defined and documentation comments are provided. 
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APPENDIX B 


HP-4l1 PROGRAM 


In the course of developing the results of Chapter III, 
it is necessary to perform a large number of simple but 
tedious calculations. the following HP-41CV program proved 
helpful. It should be easily adaptable to any small pro- 
grammable calculator. 


The program prompts the user for E/N (V-cm*) , v,/N (cm?/s) , 


N Cane ne a aoe and ae (cm/s). These data are available 
from Figure 2.2; 
The program then displays E (V/m), No (m >), a5 (m? /sm 


DTA, BETAL (laminar), BETAT (turbulent), BETAAL (fanaa 
BETAAT (turbulent) and DT in true time in microsecs. The 
minimum BETA values computed are consistent with the 
stability criterion, Equation (2.1/7). It is"q0ed  praeeiTe. 
to add about 1% to these values to ensure a smooth running 


calculation. 
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